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ABSTRACT 

In  this  paper,  the  scattered  far-field  pattern  of  a perfectly 
conducting  rectangular  cylinder  iliuninated  by  a plane  H-wave  is  obtained  by 
a modified  form  of  GTD  approximation  and  the  result  is  subsequently  improved 
in  a straightfon^ard  manner  by  eliminating  the  discontinuities  in  the  GTD 
far-field  pattern.  The  improved  scattered  far-field  pattern  is  compared  with 
results  obtained  by  other  different  approaches,  and  the  similarities  and  dif- 
ferences between  these  results  are  discussed.  Generally  speaking,  the  improved 
pattern  obtained  by  the  present  approach  is  in  good  agreement  with  the  conven- 
tional moment-method  solution.  A method  for  checking  the  accuracy 
is  also  presented  so  that  an  independent  check  on  the  satisfaction  of  the 
boundary  condition  on  the  surface  of  the  rectangular  cylinder  can  be  performed. 
The  results  of  such  an  accuracy  check  are  quite  favorable,  which  demonstrates 
that  the  improved  far-field  pattern  obtained  by  the  present  approach  is  indeed 
a close  approximation  to  the  true  solution.  The  merit  of  such  an  independent 
accuracy  checking  scheme  is  that  the  approximate  solution  obtained  can  be 
validated  without  making  comparison  with  other  methods. 
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1.  INTRODUCTION 

The  problem  of  scaccering  of  eleceromagnecic  waves  by  a rectangular 
cylinder  has  been  Invesclgaced  by  numerous  scientific  researchers  In  Che 
past;  yet  the  interest  In  this  problem  Is  still  strong,  as  evidenced  by  a 
rather  large  number  of  papers  chat  have  recently  appeared  in  the  literature 
lU.  [2],  13],  [4],  (5),  [6]  on  this  subject. 

Due  to  its  two-dimensional  nature,  the  conventional  moment  method  [2] 
can  be  applied  to  this  problem  even  when  the  cylinder  is  electrically  large, 
particularly  when  its  symmetric  properties  are  fully  exploited.  Therefore, 
it  provides  a convenient  validity  check  for  ocher  methods.  Indeed, 

Burnside  ec  al.  [1]  and  Maucz  at  al.  [3]  have  claimed  the  validity  of  their 
solutions  only  by  a comparison  of  their  results  with  the  conventional  moment 
method. 

Both  of  Che  aforementioned  papers  employ  the  so-called  uniform  geo- 
metrical theory  of  diffraction  developed  by  Kouyoumjian  et  al.  [7]  in  which 
Che  classical  wedge  diffraction  coefficient  put  forth  by  Keller  is  modified 
using  multiplicative  factors  in  such  a manner  that  one  of  these  factors  goes 
to  zero  as  the  Keller's  diffraction  coefficient  goes  to  infinity  at  one  of 
Che  shadow  or  reflection  boundaries  so  chat  the  product  remains  finite.  The 
resulting  diffraction  coefficient  is  indeed  applicable  at  all  aspect  angles. 

It  is  conceivable  that  an  additive  term  which  goes  to  infinity  as  the  Keller's 
diffraction  coefficient  goes  to  infinity,  but  with  opposite  sign,  will  work 
equally  well  as  does  the  multiplicative-factor  type  of  remedy.  Indeed,  this 
additive-term  type  of  remedy  has  been  followed  by  Lewis  et  al.  [8]  and  has 
been  investigated  by  a number  of  ocher  workers  [ sl  • 

In  this  paper,  we  Cackle  Che  scattering  problem  using  a completely 
different  technique  which  is  based  on  a representation  of  the  scattered 
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fields  la  cents  of  the  spectrum  of  the  induced  surface  current  on  the 
scaccerer  rather  chan  the  rays  emanating  from  it.  This  spectral  domain 
incerprecation  of  high-frequency  diffraction  phenomena  has  been  documented  in 
detail  in  a recent  paper  by  Mittra,  Rahmat-Samii  and  SCo  [9].  Ko  and  Miccra 
[10]  also  developed  a method  based  on  the  spectral  domain  concept  for  combining 
the  asymptotic  high-frequency  technique  and  the  integral  equation  formulation. 
The  method  has  been  applied  successfully  to  infinitesimally  chin  scattering 
objects,  e.g.,  an  infinite  perfectly  conducting  scrip  and  a finite  perfectly 
conducting  square  place  [13].  In  this  paper  we  illustrate  t.he  application  of 
this  method  to  the  problem  of  a perfectly  conducting  rectangular  cylinder. 

We  show  chat,  starting  out  with  Keller’s  wedge  diffraction  coefficient,  only 
one  iteration  gives  a far-field  pattern  chat  compares  extremely  well  with 
c.hac  obtained  from  the  moment-method  solution. 

Before  closing  this  section,  it  is  worth  mentioning  c.hac  the  CTD 
result  deviates  significantly  in  chose  directions  which  coincide  with  the 
orientations  of  the  surfaces  of  t.he  rectangular  cylinder  even  when  multiple- 
edge  interactions  are  included  in  a self-consistent  manner  by  introducing 
eight  unknown  diffraction  coefficients.  Only  when  one  further  modifies  the 
GTD  method  with  the  uniform  diffraction  coefficient  and  combines  it  with  the 
moment  method  using  a total  of  rwency-four  unknowns  will  the  result  compare 
well  with  Che  moment-method  solution  for  thirty-two  unknowns. 

2 . FORMULATION 

The  geometry  of  the  electromagnetic  scattering  problem  involving  a 
perfectly  conducting  infinite  rectangular  cylinder  of  cross  section  2a  * 2b 
illuminated  by  a uniform  plane  wave,  whose  magnetic  intensity  vector  is 
oriented  parallel  to  the  edges  of  the  cylinder,  is  depicted  in  Fig.  1.  For 
convenience  of  analysis,  an  arbitrary  incident  wave  can  always  be  decomposed 
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into  cvo  components  with  respect  to  the  z-axis,  namely,  TM^  (E-wave)  and 
TE^  (H-vave).  In  the  following  discussion  we  consider  the  H-wave  case  only; 
the  E-wave  case  can  be  solved  in  a similar  manner  by  considering  E^  ■ zEq. 

In  the  U-wave  case,  the  incident  H- field  is  given  by 

. -lk(xcos<>^+ysin«^) 

- e ^ ^ (1) 

where  the  e time  dependence  is  understood.  The  geometrical  optics 
reflected  field  can  be  derived  easily  from  (1)  once  the  geometrical  configura- 
tion is  known.  Hence,  the  geometrical  optics  part  of  the  total  field  will 
not  be  discussed  further. 

We  now  turn  to  the  diffracted  field  which  results  from  the  sharp 
edges  on  the  scattering  cylinder.  As  alluded  to  in  Section  1,  we  make  use 
of  the  Keller  wedge  diffraction  coefficient  to  find  the  diffracted  fields. 

To  make  this  discussion  as  self-contained  as  possible  and  to  introduce  the 
notations  chat  we  have  adopted  in  this  paper,  a brief  review  of  the  wedge 
diffraction  coefficient  is  in  order. 

The  geometry/  of  a perfectly  conducting  wedge  immersed  in  a uniform 
H-wave  in  the  canonical  wedge  diffraction  problem  is  shown  in  Fig.  2. 

According  to  Keller's  geometrical  theory  of  diffraction,  the  diffracted  field 
by  the  edge  of  the  wedge  can  be  constructed  by  Che  following  formula: 


H^  • edge  of  Che  wedge) 
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and  ' means  "asympcocically  equal  co."  In  the  above  formula,  nn  denoces  che 
exterior  region  while  (2  - n)^  denoces  che  wedge  angle.  Noce  chac  n is  noc 
necessarily  an  integer;  in  che  rectangular  cylinder  case,  ic  will  cake  on  che 
value  3/2.  The  subscript  w accached  co  che  angle  of  incidence  and  che 
angle  of  diffraction  > serves  as  a reminder  chac  these  angles  are  used  in 
che  canonical  problem  and  cheir  senses  are  defined  as  shown  in  Fig.  2.  This 
face  is  worth  emphasizing  since  ic  is  crucial  in  conscruccing  che  diffracted 
field  when  several  wedges  are  involved  in  a scattering  problem. 

As  shown  in  Fig.  3,  che  major  concribucions  co  che  diffracted  field 
are  from  wedges  1,  2,  and  3 since  che  chree  edges  of  these  wedges  are  being 
illuminated  by  che  incident  field  whereas  wedge  i is  in  che  dark  and,  hence, 
can  be  ignored  in  che  zeroch-order  approximation  co  che  diffracted  field. 

'<^ich  this  in  mind,  che  diffracted  far  field  can  be  wriccen  in  a concise 


6 


where  9^  - ir/4.  In  Che  above  equacion,  subscripcs  I,  2,  and  3 on  designate 
wedges  I,  2,  and  3,  respectively.  Each  of  the  from  the  individual  wedges 

is  of  the  form  given  by  (2)  and  (3)  with  proper  values  for  > and  j-  substi- 

W 

Cured  and  proper  phase  shifts  taken  care  of  due  to  the  transfer  of  the  phase 
center,  located  at  the  edge  of  the  wedge  in  each  of  the  individual  canonical 
problems,  to  the  common  phase  center,  located  at  the  origin  of  the  rectangular 
coordinate  system  shown  in  Fig.  3. 

2.1  Pole  singularicles  in  the  diffraction  coefficients 

It  is  well-known  that  the  Keller's  diffraction  coefficient  (3)  as  Intro- 
duced in  the  last  section  will  not  give  the  correct  value  for  the  diffracted 
field  at  certain  directions,  namely,  at  t.he  shadow  and  the  reflection 
boundaries  predicted  by  geometrical  optics.  One  can  easily  see  from  (3)  that 
the  denominator  of  either  one  of  the  two  terms  in  the  square  brackets  vanishes 
at  one  of  these  shadow  and  reflection  boundaries  and  gives  ficticious, 
infinite  field  values  at  chose  directions. 

Although  a number  of  uniform  theories  [7],  [8]  are  available  for  cir- 
cumventing these  difficulties,  we  suggest  a straightforward  procedure  for 
accomplishing  the  same.  We  note,  first  of  all,  chat  the  infinities  may  be 
physically  interpreted  as  chose  associated  with  physical  optics  currents  with 
semi- infinite  supports  [see  Reference  11].  Thus  we  simply  subtract  the  non- 
physical pole  singularities  and  define  a modified  diffraction  coefficient 
which  is  the  Keller  coefficient  minus  Che  pole  singularity.  With  this  modi- 
fication Che  expression  for  the  diffracted  field  from  wedge  1 (see 
Fig.  3)  becomes 
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where 


0 1 in  1 7 0 < > < -^ 


Noce  chat  Che  second  cera  in  che  square  brackets  is  che  pole  singularicy  of 
Che  first  cera  ac  che  shadow  boundary,  and  che  fourth  cera  is  che  pole  singu- 
laricy of  the  third  cera  at  che  reflection  boundary. 

Similar  expressions  for  che  diffracted  fields  H^2  con- 

structed for  wedge  2 and  wedge  3,  respectively,  in  che  appropriate  angular 
regions  (0  ^ b ^ ir  and  y £ i ^ 2it  for  wedge  2;  0 ^ b ^ tt/2  and  x ^ b i 2t 
for  wedge  3) . 

The  total  diffracted  field  may  now  be  found  by  superposing  ^jl’ 

^23  indicated  in  (4).  The  far-field  pattern  computed  by  using 
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This  plot  shows  thac 
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>g  ■ , a - b - lA,  in  (4)  is  depicted  in  Fig.  4. 

there  are  no  sporadic  variations  near  the  shadow  and  reflection  boundar  is  at 
225®,  135*,  and  315*.  This  pattern  is  already  in  very  good  agreement  with 
chat  obtained  by  using  the  moment  method,  which  is  shown  in  Fig.  7.  By  good 
agreement,  we  mean  that  all  the  locations  of  the  peaks  and  nulls  are  close 
to  the  right  positions  and  the  levels  are  more  or  less  correct.  Except  for 
the  obvious  discontinuities  in  the  pattern  at  the  angles  of  0*,  90*,  180®, 
and  270*,  which  correspond  to  the  directions  in  which  the  four  surfaces  of 
the  rectangular  cylinder  are  oriented,  it  is  really  remarkable  how  good  a 
far-cield  pattern  can  be  constructed  with  only  a straightforward  modification 
of  the  zeroth-order  GTD  solution.  A simple  physical  interpretation  for  the 
e.xistence  of  these  discontinuities  can  be  found  and  their  elimination  is 
discussed  in  the  following  sections. 

Before  closing  this  section,  it  is  .noteworthy  that,  due  to  the  com- 
plete symmetry  of  Che  rectangular  cylinder,  the  pole  singularities  in  the 
wedge  diffraction  coefficients  need  not  be  subtracted  out  explicitly  as 
done  in  (5)  before  we  use  them  to  compute  the  far  field.  In  other  words, 

Che  wedge  diffraction  coefficient  as  given  in  its  original  form  shown  in  (3) 
could  have  been  used  directly  to  compute  the  far-field  pattern.  WTiile  it 
is  admittedly  a happy  coincidence  that  the  Keller's  diffraction  coefficients 
can  be  used  directly  in  the  rectangular  cylinder  problem,  in  more  general 
cases  the  procedure  of  subtracting  pole  singularities  explicitly  from  the 
diffraction  coefficient,  as  done  in  (5),  should  be  followed  in  order  to  get 
meaningful  results  near  Che  shadow  and  the  reflection  boundaries. 

2.2  Discontinuities  in  Che  far-field  pattern  at  b “ 0,  y * ’f,  and  ^ 

In  the  previous  sections  the  problem  of  diffraction  by  a rectangular 
cylinder  was  solved  by  using  modified  Keller's  wedge  diffraction  coefficients. 


Although  the  far-field  pattern  varies  smoothly  and  remains  finite  across  the 
shadow  and  reflection  boundaries,  there  are  noticeable  discontinuities  in  the 
pattern  at  ^ » 0,  -j  , ir,  and  . These  directions  are  those  in  which  the 
surfaces  of  the  rectangular  cylinder  are  oriented.  In  the  next  section,  a 
physical  interpretation  of  the  existence  of  these  discontinuities  is  given. 
Before  doing  that,  these  discontinuities  are  explicitly  evaluated  in  this 
section. 

Let  us  refer  to  (4).  The  discontinuity  in  the  far  field  at  s « 0 is 
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Following  the  same  lines,  the  discontinuities  in  the  far  field  at  j ■ y , t, 
and  can  also  be  obtained  explicitly  [llj. 

2. 3 Physical  interpretation  of  the  existence  of  the  discontinuities  in 
the  far-field  pattern  at  b = 0,  , t,  and 


In  the  last  section,  the  discontinuities  in  the  far  field  at  d ” 0, 
and  — have  been  studied  and  documented.  The  fact  that  the  GTD  solution 
of  the  far  field  contains  discontinuities  is  well-known.  The  cause  of  the 
existence  of  these  discontinuities  is  investigated  and  a simple  physical 
interpretation  is  given  which,  in  turn,  provides  a clue  to  their  elimination. 
Note  that  in  (4)  the  far  field  is  constructed  by  using 

and  of  which  has  its  own  angular  domain  of  definition.  More 

d 3^ 

explicitly,  is  defined  on  the  angular  range  0 £ d i “j*  (see  Fig.  5a); 

outside  of  this  angular  range,  is  simply  assumed  to  be  zero.  Likewise, 
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is  defined  on  0 ^ ^ if  and  ^ ± ^ ± 2-n  (see  Fig.  5b),  and  is  on 

0 ^ ^ -j  and  If  ^ ^ 2ir  (see  Fig.  5c).  It  would  have  been  a valid  assumption 

in  the  canonical  problem  since  fields  would  have  existed  only  in  the  exterior 
region  of  the  wedge.  In  the  rectangular  cylinder  problem,  it  should  be 
realized  that  the  fields  do  not  terminate  abruptly  at  the  geometrical  planes 
represented  by  dashed  lines  in  Figs.  5a,  5b,  and  5c.  These  dashed  planes 
would  have  been  occupied  by  the  semi-infinite  wedge  surfaces  in  the  canonical 
problems.  Physically,  the  fields  should  be  continuous  across  these  dashed 
planes  in  space.  Therefore,  in  constructing  the  far  field  as  shown  in  (1), 
we  have  effectively  created  a discontinuity  at  each  of  these  dashed  planes. 

The  presence  of  these  discontinuities  is  solely  artificial  because  they  are 
created  by  our  abrupt  truncation  of  the  fields  to  the  regions  in  space 
corresponding  to  the  exterior  regions  of  the  wedges  used  in  the  canonical 
problems.  Now  that  the  cause  of  the  existence  of  the  discontinuities  in  Che 
far-field  pattern  at  <j)  = 0,  , ir,  and  has  been  found,  we  proceed  to  give  them  a 

physical  interpretation  which  in  turn  provides  a clue  to  their 
elimination. 

In  Che  following  discussion,  a typical  case  is  studied  and  the  solution 
CO  this  case  is  derived  in  a step-by-step  fashion.  The  interested  reader  is 
referred  to  Ko  and  Mittra  [11]  for  the  results  of  other  cases. 

The  typical  case  as  shown  in  Fig.  5a  is  studied.  Let  us  concentrate 
on  Che  extension  of  surface  A.  This  extension  is  represented  by  a dashed 
line  from  x=*atox  = ®aty  = b in  the  direction  iji  “ 0.  On  one  side  of 
this  plane,  there  exists  as  shown  in  Fig.  5a;  on  the  other  side  there  is 
zero  field.  By  classical  electromagnetic  theory,  we  conclude  that  there  is 
an  electric  surface  current  density  existing  in  the  dashed  plane.  Such  a 
current  radiates  in  free  space  to  support  this  discontinuity  in  the  H-field. 

In  compact  mathematical  language,  we  write 
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where  a is  the  unit  normal  vector  to  the  surface  A,  the  subscript  A refers  to 
the  extension  of  surface  A,  the  subscript  1 designates  the  wedge  under  con- 
sideration, the  superscript  d means  diffraction,  and  the  rest  should  be  self- 

,d 


explanatory.  Note  that  a factor  of  2 is  not  present  in  (7)  because  this  J 
is  simply  the  discontinuity  in  free  space  of  the  H-field  and  there  is  no 
backing  by  a perfect  conductor  in  this  situation. 

Referring  to  (5)  and  keeping  in  mind  that  the  original  diffraction 
coefficient  is  perfectly  valid  in  the  direction  j = 0,  which  implies  that 
the  procedure  of  subtracting  the  pole  singularities  need  not  concern  us 
here,  we  can  write 


xAl 


■ <» 


where 


ikp 

f^(^  = 0)  • 


o(y  - b)  , p = xe(a,») 


(8a) 


f^(<{)  » 0)  - e 


■ik  (-acosd>Q+bsin(^Q) 


.ir 


(8b) 


Note  that  in  the  above  equation,  an  approximation  of  p by  x has  been  made  and 
a constant  value  zero  for  the  variable  f has  been  assumed.  The  far-field 
expression  as  given  in  (5)  can  be  legitimately  used  to  compute  the  H-field  on 
the  dashed  plane  because  any  observation  points  in  that  plane  will  be  at  least 
at  a distance  of  2a  > 2\  away  from  the  edge  of  wedge  1. 
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The  H-field  generaced  by  is 
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ik<3  ^4 


(9) 


where  the  superscript  Ex  stands  for  "excess"  to  remind  the  reader  that  this 

field  has  been  created  artificially  by  abruptly  truncating  the  diffracted 

fields,  therefore,  it  should  not  have  existed;  hence,  the  word  "excess." 
d 
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in  (9)  is  given  by 
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By  letting  x(k  - a)  ^ J Che  integral  with  respect  to  x can  be  trans- 
formed into  an  integral  with  respect  to  t which  may  be  identified  with  the 
well-tabulated  Fresnel  integrals,  i.e.. 
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In  writing  (11),  we  have  explicitly  used  the  relation  a » ic  cos  $ and  the 
definitions  of  the  Fresnel  integrals,  i.e.. 
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Substituting  (11c)  in  (10)  and  the  resulting  expression  for  back 

into  (9),  we  have 
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where  fj^(b  = 0)  is  given  in  (8b). 

Now  we  have  the  expression  for  the  field  generated  by  as  shown  in 

(13) . The  next  step  is  to  examine  the  discontinuity  introduced  at  ^ » 0 by 
. To  this  and,  all  we  need  to  pay  attention  to  is  the  following  limit: 
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Using  (14)  and  (13),  the  discontinuity  introduced  at  4 ■ 0 by  is  round 


to  be 


r 


- 14  - 


udEx, 

h^l  (i  - 0)  - (? 


2:r) 


-ik 


i— 

ika  4 
e e 


x'Sffk!: 


-2  -[/2  - i-yl)]  ” 0)  • ^ 


i- 

4 iko 
I'  2 /TKO 


-ik(-acos®Q+bsin?Q)  (2)(-|  sin 


cos 


- cos  3 Sq 


(15) 


Ic  is  indeed  enlightening  to  see  chat  this  discontinuity  is  exactly  the 
discontinuity  given  in  (6).  The  study  of  the  typical  case  is  now  completed. 
For  the  cases  depicted  in  Figs.  5a,  5b,  and  5c,  the  interested  reader  is 
referred  to  .Ko  and  Mittra  [11]. 

Before  closing  this  section,  ic  is  worthwhile  to  recapitulate  the  main 
points  discussed.  The  discontinuities  at  ^ = 0,  ^^id  in  the  far- 

field  pattern  obtained  by  the  GTD  technique  have  been  shown  to  be  caused 
by  Che  abrupt  truncation  of  Che  diffracted  field  which  is  effectively 
equivalent  to  artificially  introducing  sheets  of  electric  surface  current 
in  free  space.  When  these  current  sheets  radiate  strongly  in  the  directions 
^ = 0,  , T,  and  , they  produce  discontinuities  corresponding  exactly  to 

chose  in  Che  GTD  far-field  pattern.  This  discovery  provides  us  with  a clue 
CO  further  improve  the  GTD  far-field  pattern.  One  simply  subtracts  out  all 
these  artificially  created  "excess”  fields  from  the  GTD  far-field  pattern  to 
eliminate  these  discontinuities.  The  improved  GTD  far-field  pattern  is 
discussed  in  the  next  section  along  with  some  comparisons  with  the  conven- 
tional moment  method  solution  and  results  obtained  by  a GTD-momenc  method 
hybrid  technique  [1]  and  a self-consistent  method  [12]. 
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3.  IMPROVED  FAR-FIELD  PATTERN  AND  COMPARISON 
WITH  RESULTS  OBTAINED  BY  OTHER  APPROACHES 

In  the  last  section,  Keller's  wedge  diffraction  coefficient  was  used 
to  construct  the  far  field  scattered  by  a rectangular  cylinder.  It  was 
shown  that  when  the  pole  singularities  contained  in  the  physical  optics 
current  with  infinite  support  on  the  semi- infinite  wedge  surface  were  sub- 
tracted from  the  wedge  diffraction  coefficient,  the  far  field  remained 
finite  and  varied  smoothly  across  the  geometrical  optics  shadow  and  reflec- 
tion boundaries.  Other  than  the  noticeable  discontinuities  in  the  directions 
in  which  the  surfaces  of  the  rectangular  cylinder  are  oriented,  that  far- 
field  pattern  was  in  fair  agreement  with  the  conventional  moment  method 
solution. 

It  was  pointed  out  that  these  discontinuities  were  caused  by  abrupt 
truncations  of  the  diffracted  field  to  the  regions  in  space  which  correspond 
to  the  exterior  regions  in  the  wedge  canonical  problems.  These  abrupt 
truncations  were  demonstrated  to  be  effectively  equivalent  to  artificial 
introductions  of  semi- infinite  electric  surface  current  sheets  in  free 
space.  An  in-depth  study  of  these  current  sheets  consequently  showed  that 
they  radiated  strongly  and  introduced  identical  discontinuities  in  those 
directions  corresponding  to  the  discontinuities  in  the  GTD  far  field.  This 
discovery  provides  a clue  which  leads  to  further  improvement  of  the  GTD  far- 
field  pattern. 

The  elimination  of  these  discontinuities  can  be  achieved  by  subtracting 

”a1  ’ ^1  ’ ”a2  ’ ^B2  ’ ^B3  ’ ^C3  ’ last  section,  from 

the  GTD  given  in  (4).  The  improved  GTD  far-field  pattern  so  obtained  is 
displayed  in  Fig.  6.  It  is  very  interesting  to  observe  that  not  only  the 
discontinuities  at  t ” 0.  and  disappear  completely,  but  also  that 

nulls,  which  were  not  present  in  the  GTD  far-field  pattern  in  Fig.  4,  are 
showing  up  in  the  neighborhood  of  (>  « it  and  . The  disappearance  of 
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the  discontinuities  is  to  be  expected  since  these  discontinuities  have  been 
examined  carefully  in  the  last  section  and  were  shown  to  be  identical  to 
those  in  the  GTD  far-field  pattern  in  Fig.  4. 

Finally,  results  obtained  by  using  (a)  the  conventional  moment  method  [2] 
with  32  unknowns,  (b)  the  hybrid  moment  method  and  GTD  technique  by  Burnside 
et  al.  [1]  with  24  unknowns,  and  (c)  the  self-consistent  method  [12]  with  8 
unknowns  are  exhibited  in  Fig.  7.  Note  that  in  Fig.  7,  the  self-consistent 
method,  which  uses  2II  higher-order  multiple-edge  interactions  in  its  formu- 
lation, still  does  not  eliminate  the  discontinuities  at  t •0°,  90*,  180®, 
and  270°.  While  the  conventional  moment  method  solution  and  the  moment 
method-GTD  hybrid  solution  deviate  only  slightly  in  the  backscattering 
direction,  the  present  method  solution  as  shown  in  Fig.  6 is  in  good  agreement 
with  them. 


4 . ACCURACY  CHECK 

In  the  previous  sections,  a method  has  been  discussed  for  obtaining  an 
improved  scattered  far-field  pattern  of  a perfectly  conducting  rectangular 
cylinder  illuminated  by  a plane  wave  whose  magnetic  intensity  vector  is 
parallel  to  the  axis  of  the  cylinder.  The  resulting  improved  scattered 
far-field  pattern  has  been  compared  with  results  obtained  by  other  approaches 
in  the  literature.  In  particular,  results  obtained  by  using  (a)  the  con- 
ventional moment  method  [2]  with  32  unknowns,  (b)  the  hybrid  moment  method 
and  GTD  technique  by  Burnside  et  al.  [1]  with  24  unknowns,  and  (c)  the 
self-consistent  method  [12]  with  8 unknowns  have  been  used  for  comparisons. 

Out  of  these  three  different  approaches,  only  the  conventional  moment 
method  has  consistently  taken  into  account  in  its  formulation  the  boundary 
condition,  requiring  the  tangential  components  of  the  total  electric  field 
to  be  vanished  on  the  surface  of  the  rectangular  cylinder.  In  Burnside's 
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hybrid  moment  method  and  GTD  technique,  the  boundary  condition  is  not  en- 
forced on  the  entire  surface  of  the  scattering  object,  as  in  the  case  of  the 
conventional  moment  method.  Instead,  point  matching  is  applied  at  the  mid- 
point of  each  of  the  pulse  current  segments,  located  near  the  diffracting 
edges  of  the  wedges  of  the  rectangular  cylinder,  and  at  some  arbitrary  points 
in  the  GTD  regions  of  the  surface  of  the  rectangular  cylinder,  in  order  to 
obtain  an  equal  number  of  equations  as  the  involved  unknowns  including  the 
current  samples  and  the  diffraction  constants.  The  final  result  is  somewhat 
dependent  on  the  locations  of  the  arbitrarily  chosen  additional  matching  points 
in  the  GTD  region  on  the  surface  of  the  scatterer.  In  the  self-consistent 
method,  no  consideration  is  given  to  the  satisfaction  of  the  boundary 
condition. 

Although  the  comparisons  of  the  scattered  far-field  pattern  with  those 
obtain<?d  by  the  aforementioned  methods  are  very  favorable,  indicating  the 
solutions  are  in  good  agreement  with  the  true  solution,  an  independent 
accuracy  check  must  be  applied  to  further  validate  the  approximate  solution, 
especially  when  there  are  no  available  results  in  the  literature  to  compare 
with.  Such  an  accuracy  check  is  often  needed  but  is  not  readily  available 
in  other  high-frequency  asymptotic  techniques. 

In  this  section,  we  present  a method  which  allows  us  to  calculate  the 
tangential  components  of  the  scattered  electric  field  on  the  surface  of  the 
scatterer  using  the  approximate  scattered  far-field  pattern.  The  accuracy 
of  the  solution  can  then  be  checked  simply  by  observing  whether  these  tan- 
gential components  of  the  scattered  electric  field  on  the  surface  of  the 
scatterer  are  equal  to  the  negative  of  the  tangential  components  of  the 
incident  electric  field  on  the  surface  of  the  scatterer.  If  the  outcome  of 
such  an  observation  is  favorable,  then  the  solution  is  good;  otherwise, 
further  improvement  is  needed. 
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4.1  Method  of  computation 


In  classical  electromagnetic  theory,  it  is  well-knovm  that  the  scattered 
magnetic  far  field  is  related  to  the  Fourier  transform  of  the  scattered 
electric  field  on  an  aperture  in  a relatively  simple  manner.  Consider  the 
aperture  shown  in  Fig.  3.  This  aperture  is  a plane  at  y ■ b,  containing 
the  surface  of  the  rectangular  cylinder  as  its  central  portion.  The 
scattered  far-field  pattern  enclosed  in  the  hemisphere  corresponding  to  this 
aperture  as  shown  in  Fig.  3 is  used  to  compute  the  Fourier  transform  of  the 
scattered  electric  field  on  the  aperture.  In  particular,  for  the  H-wave  case 
considered  in  this  paper,  the  Fourier  transform  of  the  tangential  component 
of  the  scattered  electric  aperture  field  can  be  written  as 


ik  2g(p) 


(o)  • exp(-ikb  sin  j) 


where 


g(o)  - — — ■ ■ 

>'8irkc 

Zq  is  the  free-space  characteristic  impedance 

T (?)  is  chat  portion  of  the  scattered  magnetic  far-field  pattern 
enclosed  in  the  proper  hemisphere  corresponding  to  the  aperture  under 
consideration,  and 
0 and  ? are  the  polar  coordinates. 

The  phase  factor  in  (16)  is  to  assure  chat  the  aperture  at  y = b,  but  not  at 
y * 0,  is  under  consideration.  The  tangential  scattered  electric  aperture 
field  in  the  space  domain  can  be  obtained  readily  by  an  inverse  Fourier 
transform,  viz.. 


S^(x,y  - b)  - F"^{E^[?(a)]} 


(17) 


The  tangential  scattered  electric  aperture  field  in  (17)  is  truncated  to  the 
surface  of  the  scatterer,  i.e.,  x e[-a,a],and  an  observation  can  be  made  to 
judge  whether  the  tangential  scattered  electric  field  on  the  surface  of 
the  rectangular  cylinder  is  indeed  equal  to  the  tangential  incident  electric 
field. 

The  hemisphere  shown  in  Fig.  8 is  for  the  aperture  at  y ■ b.  Similar 
hemispheres  and  apertures  corresponding  to  the  other  surfaces  of  the  rectangu- 
lar cylinder  can  be  established  and  the  tangential  components  of  the  scattered 
electric  field  on  these  surfaces  can  be  computed  in  a similar  fashion. 

4 . 2 Results  and  comments 

In  the  previous  sections,  an  improved  scattered  magnetic  far-field 
pattern  has  been  obtained  (see  Fig.  6) . The  accuracy  of  that  pattern  can  be 
verified  by  computing  the  tangential  component  of  the  scattered  electric 
field  on  the  surface  of  the  rectangular  cylinder  using  the  method  outlined  in 
Subsection  4.1.  The  resulting  tangential  component  of  the  scattered  electric 
field  on  the  surface  — x £[-a,a],  y = b — is  shown  in  Figs.  9a  and  9b. 

These  curves  were  obtained  by  using  that  portion  of  the  improved  pattern 
(Fig.  6)  enclosed  in  the  hemisphere  schematically  indicated  in  Fig.  8,  and 
by  using  the  one-dimensional  FFT  for  the  inverse  Fourier  transform  operation. 

The  magnitude  of  the  tangential  scattered  E-field  shown  in  Fig.  9a 
oscillates  around  the  constant  value  266.58,  which  is  the  magnitude  of  the 
tangential  component  of  the  incident  E-field.  More  precisely,  for  normalized 
H-wave  incidence,  |h^|  * 1,  Ie^|  » 377,  and  the  incident  angle  • ^5“; 
therefore,  the  magnitude  of  the  tangential  component  of  the  incident  E-field 
is  |E^  cos  45° I » 266.58.  The  phase  of  the  tangential  scattered  E-field  shown 
in  Fig.  9b  is  varying  linearly  across  the  surface  of  the  rectangular  cylinder 
and  can  be  readily  verified  to  be  of  a phase  difference  of  x or  180°  from  the 
phase  of  the  tangential  incident  E-field  on  the  surface.  Similar  behaviors 
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of  Che  tangential  scattered  E-field  on  ocher  surfaces  of  the  rectangular 
cylinder  can  be  obtained  and  are  not  repeated  here.  These  observations  con- 
firm that  the  improved  far-field  pattern  obtained  in  the  last  section  is 
indeed  a very  good  approximation  to  the  true  scattered  field.  Of  course, 
the  comparisons  with  results  obtained  by  other  approaches  as  done  in  the 
previous  section  not  only  further  validate  the  approximate  solution,  but 
also  demonstrate  the  effectiveness  of  the  accuracy  checking  method 
developed  in  this  section. 


5.  SUMMARY 

In  this  paper,  the  scattered  far-field  pattern  of  a perfectly  con- 
ducting rectangular  cylinder  illuminated  by  a plane  H-wave  has  been  obtained 
by  a zeroth-order  GTD  approximation  and  the  result  improved  by  a straight- 
forwardly physical  interpretation  of  the  existence  of  the  discontinuities 
in  the  zeroth-order  GTD  far-field  pattern. 

The  improved  scattered  far-field  pattern  has  been  compared  with 
results  obtained  by  other  different  approaches  and  the  similarities  and 
differences  between  these  results  have  been  discussed.  Generally  speaking 
the  improved  pattern  obtained  by  the  present  approach  is  in  good  agreement 
with  the  conventional  moment-method  solution  with  32  unknowns. 

An  accuracy  checking  method  has  also  been  presented  so  that  an  indepen- 
dent check  on  the  satisfaction  of  the  boundary  condition  on  the  surface  of 
the  rectangular  cylinder  can  be  performed.  The  results  of  such  an  accuracy 
check  are  quite  favorable,  which  demonstrates  that  the  improved  far-field 
pattern  obtained  by  the  present  approach  is  indeed  a close  approximation  to 
the  true  solution.  The  merit  of  such  an  independent  accuracy  checking  scheme 
is  that  the  approximate  solution  obtained  can  be  validated  without  making 
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comparisons  with  other  methods.  Although  not  needed  here,  the  spectral 
approach  to  diffraction  can  be  employed  in  conjunction  with  the  integral 
equation  for  the  surface  current  to  systematically  improve  the  solution 
where  such  an  improvement  is  needed. 
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OBSERVATION  POINT 


Fig.  2.  Geometry  of  a perfectly  conducting  wedge  immersed  in  a 
uniform  H— wave  in  the  canonical  wedge  diffraction 
problem. 
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Fig.  5a.  contribute  to  the  discontinuities  in  the  far 
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field  at  # ’ , cfi  = 0,  respectively. 
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Fig.  7. 


db  plot  of  scattered  far-field  patterns  for  the  rectangular 
cylinder  calculated  by  using  the  moment  method,  GTD  and  hybrid 
technique;  |i  a“b*lX. 
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